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Abstract 

We study the energy distribution of <j> in pure penguin induced B — > X 3 <p 
taking into account the fermi motion of b inside B meson for b — > s(j) and also 
modification due to gluon bremsstrahlung process b — > s<pg. We find that the 
contribution to B — ► X s (p from b — > scftg is less than 3%. This study provides 
a criterion for including most of the 0's produced in a penguin process. 
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Rare B decays, particularly pure penguin decays, have been a subject of considerable 
theoretical and experimental interest recently [1] . The photonic penguin processes have been 
observed by the CLEO collaboration || in both the exclusive mode B — ► K* / y and in the 
inclusive mode B — > X s j. The Standard Model (SM) is consistent with experimental data 
||. A signature of pure penguin hadronic processes are the exclusive modes B — > K(f>, K*(p 
@] etc. or the inclusive mode B — > X s cf) |5||| and other modes resulting from processes like 
b — > sss. The search for exclusive processes has not yet led to a definite observation. The 
inclusive mode with a larger branching ratio would be a complementary way of searching 
for penguin processes. At the quark level B — ► X s <f) results from b — > s<p, just as B — > Xg^ 
results from b — > sj. In both cases the energy spectrum of or 7 are not monenergetic as a 
result of two effects. First, because b quark is in the B meson, its fermi momentum smears 
the energy spectrum of <fi or 7. Therefore the distribution of energy depends to some extent 
on the choice of the wave function. Second effect arises from the process b — > sg'y in the 
photonic case, and b — > sg<p in the hadronic case. The photonic case has been discussed in 
a series of papers by Ali and Greub [7j]. They find that the dominant contribution to the 7 
spectrum comes from the wave function effect. We shall perform a similar calculation for the 
hadronic case. The wave function effect is treated with a Monte-Carlo simulation of decays. 
The gluonic correction is carried out in a simple effective Hamiltonian approximation. We 
find that the second effect is negligible in our case. 

The QCD corrected H^ B= i relevant to us can be written as follows ||: 

#ab=i = ^[K&KUeiO? + c 2 o u 2 ) + v cb v: s ( Cl oi + c 2 o c 2 ) - v tb v* £ aOi] + h.c. , (1) 

where the Wilson coefficients (WCs) q are defined at the scale of /1 ~ m h \ and Oj are defined 

as 

Of = s a 7 M (l - 75)<?/3g/37 M (l - l5)b a , 0\ = sj^l - 75)gg7 M (l - lb)b , 
3 ,5 = *7m(1 - 7s)& ^7m(! =F 7s)</ , Qifi = s a 7 M (l - 75)6/3 v'fsl^i 1 =F 7sK , 

q> q> 

3 3 
07,9 = g^^ 1 -75)bJ2e q ,qY{l ± 7sV , <5s,io = g s a 7^(l - 75)^^ e g'^7 M (l ± 7s) 9^ • (2) 
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The Wilson coefficients at /i = mj at the next-to-leading order have been evaluated in 
Refs. @,|,|. For m t = 176 GeV and a s (m z ) = 0.117, we find 

Cl = -0.307 , c 2 = 1.147 , c 3 = 0.017 , c 4 = -0.037 , c 5 = 0.010 

c 6 = -0.045 , c 7 = 1.2 x 10~ 5 , c 8 = 3.8 x 10~ 4 , c 9 = -0.010 , c 10 = 2.1 x 10~ 3 . (3) 

We shall consider the effect due to fermi momentum in b — > s</> process, which we assume 
has the same energy distribution as 5 — > X s 0, in section (a) below. In section (b) we 
shall consider the effect due to b — > s0g process. 



a. b —> S(j) 

Using Hab=i in Eq.(l), we obtain the decay amplitude for B — > X s (p 

A(b - S 0) = -as 7M (l - 75 )ty" , (4) 

where e M is the polarization of the <fi particle; a = (g^GpVtbV^/ v2) [C3 + C4 + C5 + £(03 + C4 + 
c 6) ~ ( c 7 + c 9 + c io + £( c s + c 9 + c io))/2] with £ = 1/N C , and iV c is the number of colors. 
The coupling constant is defined by < 1 s'y^s 1 >= ig^. From the experimental value 



for Br((j) — > e + e ) \TU\, we obtain g\ = 0.0586 GeV . The branching ratio for b — ► s0 is 
predicted to be 1.7 x 10" 4 @ for a s (m z ) = 0.117. 
The decay rate is given by 

r ( 6^) = ^A^[i + f4(i-4 + 4)], (5) 

where \j = (1 — m 2 Jm 2 — rrif/ml) 2 — 4m 2 m 2 /mf. 



To study the energy distribution of 0, we adopt the model in Ref. |TT| in which the 
b quark is not at rest inside B but with a fermi momentum p& according to a Gaussian 
distribution, 
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where the parameter pf is determined from experimental data to be between 0.21 to 0.39 
GeV ||12|| . The b quark mass expressed in terms of the B meson mass and the spectator 
quark mass m q in the rest frame of B, is given by 



m t = m 2 B + m 2 -2m B ^p 2 b + m 2 q . (7) 

In the rest frame of the B, the b — > s<p decay width T{m b ) is given by (m b / E b )T(b — > s<p). 
Its contribution to the decay width T S( p for B — > X s (p is averaged over all allowed momenta 
Pb- We have 



r s<f> (B -> X s( j)) = [ PmaX $(p b )f b T(m b )d\p b 
J o 



P„ 



(m 2 B + m - (m + m s ) 2 ) 2 

\ a 2 m a ■ ( 8 ) 

Due to the finite momentum distribution, the energy of <fi from b quark decay is no longer 
monoenergetic, instead there will be a distribution. The 4> energy spectrum generated from 
a Monte-Carlo simulation of the decay is shown in Figure 1. In the figures we have used the 
constituent mass of 0.3 GeV and 0.5 GeV for the spectator quark and the s-quark, respec- 
tively. From Fig.l, we see indeed that there is a spread in the <ft energy with the maximum 
located at about 2.55 GeV. 



b. b — > sc 

The quark level effective Hamiltonian responsible for b —>■ s<pg is complicated. We use the 
simplified effective Hamiltonian in Eq. (4) to obtain the <fi energy distribution for the process 
b — > s<pg by attaching a gluon on either the initial b or the final s quarks. This is expected 
to be a good approximation because the dominant effect comes from the bremsstrahlung of 
gluon emission from the external light quark. We obtain 

,s llM {2pl- ft g Y)T a (l- ls )b 



A(b -> s(pg) = ag s (- 



{Pb ~ P g ) 2 - 



mi 



{Ps + P g ) 2 -m 2 s v 
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where pb ,p s , and p g are the b-quark, s-quark and gluon momenta, respectively. From above 
we have the following (f> energy spectrum 

dT(b^>s4>g) \a\ 2 a s N 2 - 1 f*™* , , 1 v . 2 2 1 + 



/ ^Try (4 + (1 + F) +(!-^ 

" tm.in. -L I ■* 



2 '1 - 2/4 - * - Y + XY - 2 f S s \^-)( (1 ~f s)2 -2^ + 1 + /i 2 )] , (10) 



' X 2 (1 + Y) K ^ ^ S 1 + Y JK fi 2 

where 

s + t — m 2 — rr& s + t — m 2 — m\ 



X= == v - , Y 

s = ml + ml- 2m b E (j) , (x s ^ 



'2 ' i 9 i 2 ' 

s — t — + m 

"7 S ,<£ 
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The energy distribution in Eq.(10) has the well-known infrared divergence due to the zero 
mass of the gluon. To regulate the infrared divergence, we assign an effective gluon mass of 
about 2777.^ which represents the lowest invariant mass of the gluon. The <fi energy distribution 
for b — ► s<pg is shown in Figure 2. Here we have neglected the effect due to non-zero pb 
discussed in the previous section which is small and approximated the b —>■ S(pg contribution 
to B — > X s cf) by Eq.(10). We find that the effect of b — > stfig on B — > X s <f) is small because 
BR{b — > s(j)g)l 'BR(b — > s<fi) is only about 3%. The total energy distribution is shown in 
Figure 3. 

The spectrum of <fr should approximate the spectrum that arises from the decays B — > 
Kn(f), K7nr<f), etc. Of course the monoenergetic 0's that arise from two body modes like 
B — > K(j) or K*<p are included in an average sense. The specific two body modes are not 
expected to be more than 10% of the inclusive X s (p production ||. In this paper we have 
not discussed the energy spectrum from the decay of the dominant non-penguin processes 
which are expected to have a much softer spectrum since they always arise from decay 
of charmed states. We assume that this experimentally well known contribution has been 
substracted in the region of interest. If in addition to the selection criterion on discussed 
in this letter, X s will be experimentally shown to include an odd number of kaons, then the 
penguin process will be even more enhanced Q. 
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FIG. 1. E$ distribution for b -> s<j> for P/ = 0.3 GeV. 
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FIG. 2. £^ distribution for b — > s^. 
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FIG. 3. E$ distribution for b — > s0 + s0g. The solid, doted and dashed lines are for pj= 0.39 
GeV, 0.30 GeV and 0.21 GeV, respectively. 
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